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ABSTRACT 

Aims. Already from the initial discoveries of extrasolar planets it was apparent that their population and environments are far more 
diverse than initially postulated. Discovering circumbinary (CB) planets will have many implications, and in this context it will again 
substantially diversify the environments that produce and sustain planets. We search for transiting CB planets around eclipsing binaries 
(EBs). 

Methods. CB-BLS is a recently-introduced algorithm for the detection of transiting CB planets around EBs. We describe progress 
in search sensitivity, generality and capability of CB-BLS, and detection tests of CB-BLS on simulated data. We also describe an 
analytical approach for the determination of CB-BLS detection limits, and a method for the correct detrending of intrinsically-variable 
stars. 

Results. We present some blind- tests with simulated planets injected to real CoRoT data. The presented upgrades to CB-BLS allowed 
it to detect all the blind tests successfully, and these detections were in line with the detection limits analysis. We also correctly detrend 
bright eclipsing binaries from observations by the TrES planet search, and present some of the first results of applying CB-BLS to 
multiple real light curves from a wide-field survey. 

Conclusions. CB-BLS is now mature enough for its application to real data, and the presented processing scheme will serve as the 
template for our future applications of CB-BLS to data from wide-field surveys such as CoRoT. Being able to put constraints even 
on non-detection will help to determine the correct frequency of CB planets, contributing to the understanding of planet formation 
in general. Still, searching for transiting CB planets is still a learning experience, similarly to the state of transiting planets around 
single stars only a few years ago. The recent rapid progress in this front, coupled with the exquisite quality of space-based photometry, 
allows to realistically expect that if transiting CB planets exist - then they will soon be found. 

Key words, methods: data analysis - stars: variables: general - stars: planetary systems - occultations - binaries: eclipsing 



1. Introduction 

Bound binary stars are one of the most common environments in 
the Galaxy. Studying planet formation and evolution without ac- 
counting for binary (and multiple) star systems is surely incom- 
plete - so we aim to try to fill-in observational data on these types 
of systems. To help achieve this goal, a Binaries Working Group 
was recently formed within CoRoT Exoplanets Science Team of 
the CoRoT space mission (Baglin et al 2006) to coordinate the 
different searches for planets in binary systems, and we report on 
one of the activities within the Working Group, namely - transit 
searches using the CB-BLS algorithm (Ofir 2008, hereafter pa- 
per I). These efforts are now all the more relevant following the 
recent detection by Lee et al (2009) of an eclipse- timing signal 
compatible with two circumbinary planets around HW Virginis. 

Planets around eclipsing binaries may be expected in signif- 
icant numbers, since CoRoT finds close and deeply eclipsing bi- 
naries in about 1% of the stars surveyed (Almenara et al 2009), 
leading to a sample of about 300 close binaries from its first year 
of operation. Assuming that a planetary orbit would be roughly 

* Based on observations obtained with CoRoT, a space project op- 
erated by the French Space Agency, CNES, with participation of the 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, 
Germany and Spain. 



aligned with a binary component's orbital plane, EBs constitute 
a sample with planetary orbits being preferentially aligned to 
display transits across the host stars. 

There are many detection techniques for planets in bina- 
ries: Some are common to planets around single stars (such 
as: radial velocity (RV), transits, astrometry and microlensing 
[Muterspaugh et al 2007, Lee et al 2008]) and some are unique 
to binaries (such as: eclipse timing [Deeg et al 2000, Lee et al 
2009] and gravitational waves [Seto 2008]). Generally speak- 
ing, in each geometry (i.e., circumprimary or circumbinary or- 
bits) techniques that are in principle the same take on different 
emphasizes. 

In this paper we will focus on transits of CB planets, 
and specifically transiting CB planets around eclipsing binaries 
(EBs). The transit signal from CB planets is not periodic, hence 
common techniques searching such signals are of low efficiency 
and special search algorithms have to be employed, such as CB- 
BLS [paper I] or TDA [Doyle et al 2000]. In the following we 
list the accumulated additions to CB-BLS since paper I in §2, 
describe blind tests to CB-BLS in §3, and show examples from 
the first application of CB-BLS to real data in §4, and conclude. 
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2. Accumulated Additions to CB-BLS 

The following describe additions to the CB-BLS implementation 
starting from the version used to generate the results of paper I 
(internally designated as version 0.51) till now (version 0.83). 
While the basic idea of orbit fitting remains unchanged, the spe- 
cific implementation has improved significantly in search sensi- 
tivity, generality, speed and results analysis. To avoid extensive 
repetitions, we assume in the following that the reader is familiar 
with paper I. We note that we plan to release the CB-BLS source 
code (written in MATLAB) in the future. CB-BLS now: 

- Allow the inclusion of the surface brightness ra tio ( /) in the 
input, which increases sensitivity. See more in |2.1| 

- Allows natural use of the more accurate Roche-lobe geome- 
try. All the user has to do is to specify the surface potentials 
£l\ 2 - which are the direct outputs of EB modeling - in- 
stead of the radii. CB-BLS then computes the 3D shape of 
the surface, rotates the shape to the binary inclination, and 
computes the silhouette of the sky-projected shape at each 
binary phase. CB-BLS allows single- surface binaries (con- 
tact binaries) too - which are common. We note that this 
feature was forecasted in paper 1 (§3.2) and is implemented 
along these lines. We comment that using this feature slows 
CB-BLS significantly. 

- Accounts for EB inclination. The EB orbital inclination is 
calculated anyhow during the EB modeling, so the 2D (in- 
stead of ID, as in paper I) sky-projected orbit for each com- 
ponent is computed - and the 2D distance from the test plan- 
etary model found. The planetary orbit is still assumed edge- 
on, but this is a good approximation since the planet is farther 
away from the baricenter than the EB components, and so it's 
alignment requirements are more stringent to begin with. 

- Enables the simultaneous solution of multiple light curves of 
the same system, each with it's own EB model - which also 
allows one to use multi-band data in a single run. Since after 
regularization (see paper I §3.3) the transit signal should be 
achromatic - the CB-BLS statistic can be calculated using 
all the regularized residuals - regardless of band. We note 
that performing a multi-band solution and using the (single) 
surface brightness ratio / are mutually exclusive. 

- Includes the directional correction (Tingley 2003) - ignores 
negative depths. 

- I nclu des a set of utilities: detection limits estimator (see 
§ |2.2| ), CB transit predictor, non-linear optimization. 

- Has better visualization of results for analysis, and has im- 
proved robustness, better speed optimization, more input 
control, bug fixes etc. 

- Allows one to output all the computed CB-BLS values (and 
not just the maximal one for each orbital period) to properly 
account for the sample of ?tted models [Ofir 2009]. 

We stress that CB-BLS is a detection algorithm, and as such 
- concisely avoids many complicating elements that are known 
to affect the shape of the transit light curve, including: stellar 
limb darkening, motion of EB components along the line of 
sight (the very basis of RV and eclipse timing techniques), plane- 
tary orbital eccentricity and obliquity, non-Keplerian (especially 
Newtonian) motions, etc. The reason none of these is accounted 
for is discussed in the last paragraph of the following section. 



2.1. CB-BLS with a single depth 

If the EB surface brightness ratio / is known (and / is a natural 
output of most, if not all, EB modeling tools) then the depth 



ratio between transits in front of the primary and the secondary 
EB components will also be /, and eq. 1 in paper I can then be 
re- written with just a single fitted depth as (note the last term): 

d = Yj w ^ xk - h ?+Tj Wi(Xi - L)1 + Z w j (x j ~ JL)1 (1) 



Here we approximate L 2 from paper I to JL, which is an ap- 
proximation to the correct value of J(L + H). The approximation 
is good as long as \H\ <sc \L\, which is the case for low duty-cycle 
signals such as transits. As in BLS (Kovacs, Zucker and Mazeh 
2002, hereafter KZM) and CB-BLS paper I, minimization allows 
to analytically compute H (unchanged) and L (revised): 



H 



-(si + Sj) 
1 - (n + r 2 ) 



+ JS2 

V\ + J 2 r 2 



(2) 



where the ri ?2 and s\^ 2 are simple sums of the data and the 
weights as in paper I. Plugging these to D gives: 



d = Z Wnx » ~ i 



(si + s 2 ) 2 (si + Js 2 ) 2 



■ (n + r 2 ) n + J 2 r 2 



(3) 



Similarly to BLS and to paper I, the first term on the right 
hand side is constant, leaving the rest as SR - the new CB-BLS 
statistic (times -1): 



SR 



(Si + s 2 f (si + Js 2 f 



1 - (n + r 2 ) n + J 2 r 2 



(4) 



The use of /, together with the use of Roche geometry, EB 
inclination and the directional correction, are all attempts to use 
all available knowledge about the host EB and the expected tran- 
sit signal shape to reduce the number of free parameters. For this 
reason CB-BLS uses the simplest possible box-like transit model 
and concisely avoids even slightly more complicated models that 
requires an additional free parameter (such as the trapezoidal 
transit model, limb darkening, or the other effects mentioned 
at the end of the previous section). Note however, that the EB 
model itself can be as complicated as needed - as long as it stays 
a pre-processing step to CB-BLS it does not add free parameters 
to the CB-BLS fitting process. 

2.2. Estimating CB-BLS detection limits 

Giving robust estimates for detection limits of CB planets is dif- 
ficult. A full analysis needs to account for all the effects known 
from planets around single stars, and for a number of additional 
significant sources of noise, for example: EB modeling errors, 
which are system-dependant and introduce systematic errors of 
unknown magnitude. We believe that such a full analysis can 
only be done with extensive simulations on a system-by-system 
basis. Still, as long as the approximations assumed in Paper I still 
hold, especially that the CB planet moves in a circular and edge- 
on orbit and that the noise is purely Gaussian, one can formulate 
an analytic detection limit. Fig. [T] shows the side- and front- (or 
observers-) views of the geometry of one component of an EB, 
explaining the different symbols used below. 

We begin by considering the effective transit Signal/Noise 
ratio for single stars defined by KZM: a sing i e = ~ ^[Nd where S 
is the transit depth, cr is the dispersion of the LC, N is the num- 
ber of data points and d is the transit duty cycle - the fraction of 
in-transit data points. If the surface brightness ratio / is known 
then one can extended a to transiting CB planets; we remind the 
reader that after regularization the transit depths are only two 
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discrete values with a ratio of /. We therefore rearrange to get 
5 ^/Nd = acr. In order to account for the fact that two stars can be 
transited, we change the left side to 6 ^N(d\ + d 2 ) where d\ , d 2 
are the individual on-transit duty cycles of the components, and 
in order to account for the different depths of the two compo- 
nents, we make the following changes: JN(d\6^ + to give 
a final definition of 



Secondary 



a CB = ~ jN(di + d 2 J 2 ) 



(5) 



where S is the regularized depth of the primary transit. When 
investigating detection limits one must first define what "de- 
tection" is. Using the general guidance of Fig. 6 on KZM for 
& single, corrobated by our experience from the Blind Tests of $3] 
(where e.g. 'Test3' with a = 11 was near the limit of detectabn- 
ity, see Fig. [4]), we estimate that robust detections usually have 
a > 10-12. We then use the acB definition to get 



6 = 



cra C B 



yfmli + d 2 J 2 ) 



(6) 



We now need to find d\, d 2 . For a particular system the duty cy- 
cles depends significantly on the exact relation between the or- 
bits of the two binary components and the CB planet - a small 
change in the ephemeris can lead to the drastically longer and/or 
shorter transits. We therefore wish to avoid this problem and 
consider only the average detection power at a particular plane- 
tary orbital period. In analogy to the transit duration calculation 
around single stars, the average duty cycle for each of the two 
components is d - where L is the transit chord lenght of a 
planet whose orbital plane is spanned by the line of sight (e.g. i p 
= 90°) and the line of the component's nodes, as shown in Fig. 

m 

Since in Fig. [T] the height h can be found from both the side 
view as: cos (4) = \ and from the front view as: sin(/?) = | one 
gets the relation: 



L = 2Rcos(fl) = 2RJl- 



acos(ib) 
R 



(7) 



where a is the radial projection of the component's distance 
from the center-of-mass (for the definition of a, ft, h see Fig.[T]). 
Note that L can vanish completely when acosfe) > R (i.e., the 
planet passes above/below the binary component without transit- 
ing it). We note that a should scale differently for the primary and 
secondary components according to their mass ratio q, and that a 
is itself a function of binary phase: each component will have the 
maximum chord length of 2R at the orbital nodes (where a = 0), 
and lower values - as low as L of Eq. [7]- at other times. To bet- 
ter estimate L the time-average a for each component should be 
determined, which for circular orbits (and L > at all times) is 

_ 9 a ( P \2/3 

a = -a/,. Using Kepler's third law ^ = y^J and substituting 
back one finds that: 



d = 



R_lPb 
na b \P P 



2/3 



l- 



a cos(4) 
R 



(8) 



For each of the binary components the R and a will be different, 
giving d\ 2 for the two components - to be substituted back into 
Eq.0 

Since i, J, Pb,N,cr and the ratios: R/at and R/a are all given 
by the EB model and the LC, we now have a detection limit S 




Front view 




/ ■ Primary 
Center of mass x : 

Plane of the sky 

Fig. 1. Geometrical layout of one component (marked 
"Secondary") of an EB for the estimation of CB-BLS detection 
limits The primary is shown only for easyer understanding of 
the layout. Symbols used are z& for the binary orbital inclination, 
a for the radial projection of the component's distance from 
the center-of-mass, h for the maximum projected distance of 
the component above the line sight (or more exactly: above the 
plane containing both the line of sight and the line of nodes), 
and the angle ft between a stellar radius and the chord L across 
the stellar disk cut by the above plane (boldfaced). 



as a function of only our empirical choice of acB and the EB 
geometry model - depending only on the planetary period P p . 
We note that the derived detection limit is for the best-detectable 
configuraton of planetary orbits, which are those that are copla- 
nar with the line of sight and at the same position angle as the 
EB's orbit. Finally, we remark that fulfilling the Gaussian noise 
requirement may be achieved by quatifiying the correlated noise 
(e.g. as cr r in Pont et al. 2006) and adding it to the cr used in the 
formulation above. 



3. Blind tests to CB-BLS 

CB-BLS was subsequently blind-tested using a detached EB ob- 
served during the initial run of CoRoT (CoRoT ID 102806577). 
Note that the data were not detrended in any way. Notably, 
CoRoT's orbital period was not removed from the data. 
Therefore, the tests are done in the presence of strong systematic 
noises. The "examiner" (HD) created the test data in four steps: 
(1) Removal of a constant slope apparent in the data and conver- 
sion to relative flux units (F - Fq)/Fq, where Fq is an average of 
the out-of eclipse flux. The test LCs contain 7213 data points. (2) 
Generation of a model light curve (LC) of the EB and removal of 
this model from the data. (3) On the residuals, suppression of the 
nois^]in a band-pass of periods between 5 h and 12 h. Random 
noise of similar amplitude was then inserted in that band-pass. 
This band-pass corresponds to periods a few times longer than 
the typical duration of transits events; their detectability should 
therefore not be affected. Also, periodograms of the lightcurves 
before and after this operation were compared in order to assure 
that the overall noise characteristics did not change significantly. 
(4) Adding to the previous data a model LC that consists of the 
EB model of step (2) together with a simulated circumbinary 
planet, using the UTM transit simulator [Deeg 2009]. Steps (3) 
and (4) were repeated six times, simulating different planets in 
step (4). 

The CB-BLS "solver" (AO) then received the six test LCs 
designated TestO through Test5, and parameters of the physical 
model of the system. The solver then independently solved the 
LC using JKTEBOP [Southworth et al 2004a,b] and applied 



1 This step disables 'detections' based on simple subtractions be- 
tween the various test light curves, all of which had been generated 
from the same original. For test-sets based on different originals, this 
step would not be necessary 
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Table 1. Main EB- and CB planet- parameters that affect the 
planet's transit. 



Parameter 


Value 


Unit 


Period 


3.6670288245 


[d] 


Inclination 


86.6699982 


[deg] 


Ri 


1.856 


[Re] 


R 2 


0.6838 


[Re] 


U 


14.620 


[L Q ] 


L 2 


0.2500 


[As,] 


RMS of residuals 


~ 1.3 


[mmag] 


TestO R p 


1.09 


[Rj] 


Testl R p 


0.95 


[Rj] 


Test2 R p 


0.78 


[RA 


Test3 R p 


0.79 


[Rj] 


Test4 R p 


0.95 


[Rj] 


Test5 R p 


none 


[Rj] 



CB-BLS. For the CB-BLS analysis, only the mass ratio was used 
out of the full physical model. As written in paper I, if the mass 
ratio is not known - it can simply be searched-on as a free pa- 
rameter. Figure [2] depicts the EB itself with the JKTEBOP model 
(from the Test5 LC), and also the six residuals LCs. Note that the 
residuals have some structure, meaning, they have some corre- 
lated noise - as expected in real data. In Table [T] we list some of 
the EB's and simulated planets' physical model parameters that 
most affect transit detection. Note that the RMS of the residuals 
is about 1 .3mmag. We comment that each of these CB-BLS anal- 
yses take approximately 2min on a 3 GHz PC (using the natural 
parameter resolutions presented in paper I), so processing time 
is neglible, given the small number of targets relative to constant 
stars. 

The six tests signals were such that TestO was relatively 
easy and not blind (it was used primarily to make sure the test 
procedures are working) - but still the simulated system was 
a physically -possible system, while Tests 1-5 were completely 
blind. Fig. [3] shows the resultant CB-BLS periodograms and 
Table [2] gives the resultant best CB-BLS fits values. All six test 
LCs were correctly identified, where notably Test5 was correctly 
identified only when the single-depth SR of jj |2.1| was used. Note 
that the simulated planets had their the inclination and position 
angle varied to some extent and that their LCs also account for 
limb darkening - but since these none of these features are deriv- 
able from CB-BLS they are not given in this comparison. 

While simple phase-folding the LC will not separate well in- 
and out- of transit points (which is the motivation behind CB- 
BLS), a "folded" graphical representation of the CB planet's 
transits is possible (and instrumental for signal analysis): in- 
transit points are expected when the sky-projected distance be- 
tween the planet and one of the components is closer than 
that component's radius. For that reason the data can be folded 
against the minimum distance of the planet from each of the EB 
components - and that distance is scaled by that (i.e., the closer) 
star's radius. The result is that all in-transit points are expected 
below scaled-distance of 1, and all out-of-transit points are ex- 
pected at scaled-distances grater than I. Fig. [5] shows the LC of 
Test3 folded in this way, which can also serve as an illustration 
for the ability of CB-BLS to detect even shallow transits signifi- 
cantly. Note that this simulated planet caused only three primary 
transit events, where a primary transit is only ~ 1.4 mmag deep 
(to be compared with RMS dispersion of 1.3 mmag). This depth 
corresponds to a planetary radius of 0.737?/ in an EB with a pri- 
mary component with a radius of 1.8567?©. Scaling this depth to 
\R Q (and ignoring the much-fainter secondary component) gives 



, Test LC Detection Limits For Different ot, 

x10" 3 C 

1.8 r 




"5 10 15 20 25 30 

Planetary period (d) 

Fig. 4. Detection limits for a = 10, 1 1, 12 of the CoRoT test LC 
(sold lines) and the Test3 signal as it was detected (from Table 

0. 



0.397? / - smaller than all known transiting planets but the unique 
super-Earth CoRoT-7b [Rouan et al , 2009]. As seen on Fig. 
[3] the Test3 LC is significantly detected, but is not far from the 
limit of reliable detecti on. I ndeed, one can also apply the detec- 
tion limits procedure ( §2.2\ to the CoRoT test LC too - and Fig. 
[4] shows these for various acB values with the detected signal for 
reference - and the detected signal is indeed just above or just 
below the detection limit, depending on the choice of a cb- With 
the successful detection on all test signals and the agreement be- 
tween the theoretical detection limits and actual detection, we 
conclude that CB-BLS passed the blind tests successfully and is 
ready for use on real data. 



4. Application of CB-BLS to real data 

4.1. TrESLyrl field 

The photometry of the TrES survey Lyrl field [O' Donovan et 
al 2006, Dunham et al 2004] is freely available at the NASA 
NStED website^] but this photometry is after TrES's detrending, 
and so can't be used as-is for our purpose. We therefore asked 
the above authors for the raw photometry of the same data, and 
we process below only the raw Sleuth data - which contributes 
~ 11,000 data points of the total ~ 15,000 of the TrES Lyrl 
dataset. 

Here we present our first application of CB-BLS to pho- 
tometry from a real wide- field survey and the following pre- 
processing scheme is the template for our future applications of 
CB-BLS. The main difference between the scheme below and 
the corresponding pre-processing steps of other transit searches 
is the ability to not only detrend intrinsically-constant stars but 
also to detrend variable stars by iteratively detrending their resid- 
uals. This approach was first proposed by Kovacs et al (2005) 
for their reconstructive TFA algorithm of transit signals, but here 
we aim to de-trend EBs, which show large variety of variations, 
and so - the bin-averaging method proposed by Kovacs et al 
(2005) is ill-suited for the current problem. The pre-processing 
has four main stages: 



2 http://nsted.ipac.caltech.edu/ 
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Fig. 2. Panel A: LC of the EB used in the blind tests (dots) and its JKTEBOP model (line). This LC is the Test5 LC which has no 
planet added to it. Panel B: the residuals of the six test LCs from TestO (top) to Test5 (bottom). Each successive residuals LC was 
shifted by 1% to aid visibility. 



Table 2. The input values of the simulated CB planets ("Examiner" rows) and their respective solution ("Solver" rows). The possible 
detection levels are: "+++" (Strong and secure), "++" (Probable. Noisy), "+" (Weak detection) and "no planet". The depth of the 
secondary transit is omitted since it is simply / times the primary transit. 



Test 


Attribute 


Detection 


Period 


Baricenter 


Depth 


In-transit 


& single 


Comment 


dataset 






[d] 


crossing 
[CoRoT JD] 


(primary) 
[mmag] 


points 




TestO 


Examiner 
Solver 


+++ 


13.54 
13.545 


2594.54 
2594.525 


4.683 
4.768 


148 


44.0 


Not blind. 


Testl 


Examiner 
Solver 


+++ 


12.673 
12.667 


2563.567 
2595.251 


2.749 
2.494 


145 


22.8 


A bug (since then fixed) caused 
the detected reference time to 
be when the planet crossed the 
baricenter behind the EB system, 
causing a half-integer difference 
(2.5 periods) from the correct one. 


Test2 


Examiner 
Solver 


+++ 


13.08 
13.091 


2594.54 
2594.513 


1.853 
1.560 


248 


18.6 




Test3 


Examiner 
Solver 


+++ 


17.748 
17.738 


2597.54 
2597.571 


1.901 
1.349 


132 


11.7 




Test4 


Examiner 
Solver 


+++ 


12.673 
12.652 


2563.567 
2595.295 


2.749 
2.579 


143 


23.4 


See comment to Testl 


Test5 


Examiner 
Solver 


no planet 
no planet 














A very weak signal was ini- 
tially suspected before using the 
single-depth analysis. That pe- 
riodogram peak completely dis- 
appeared when this feature was 
added and used. 



1. Determination of the systematic effects: This is very 
similar to what is done routinely on transit surveys. Firstly, out- 
lier data points are removed form all light curves by sigma clip- 
ping each light curve around a small-window (5 points) median 
filter. Next, we wish to determine a few S YSREM effects [Tamuz 
et al 2005] - but we must make sure that variable stars are not 
part of the set of stars that is used to determine the effects. To 
that end we use the Alarm statistic [Tamuz et al 2006] as a gen- 
eral variability statistic (used similarly to the Stetson J statistic 
[Stetson 1996]) in the following manner (Note that in general 
Alarm will have a smaller value for less- systematic input): 1) 



We compute Alarm for all stars 2) we determine the sub- set of 
constant stars using an Alarm maximal value that corresponds 
to the bulk of the stars (i.e., many stars that change similarly, so 
they probably don't have intrinsic variability). See Fig. [6] for an 
example. 3) The SYSREM effect(s) are computed using only this 
subset of constant stars (one effect in the above example) and the 
correction is applied to all the stars, 4) Alarm is re-computed. 
Since now the bulk Alarm distribution is narrower (e.g., main 
panel of Fig [6]) it allows to better filter true variables - so us- 
ing a lower Alarm maximal value the set of constant stars is re- 
determined. Steps 2-4 can be repeated several times till the set 
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SDE of TestO 



SDE of Testl 



SDE of Test2 




0.04 0.06 0.08 0.1 0.12 



0.04 0.06 0.08 0.1 0.12 0.14 



SDE of Test3 



SDE of Test4 



SDE of Tests 




0.04 0.06 0.08 0.1 0.12 0.14 



0.04 0.06 0.08 0.1 0.12 0.14 
Orbital frequency 



0.04 0.06 0.08 0.1 0.12 0.14 



Fig. 3. The straightened periodograms of the blind tests — to be read together with the input and output values specified on Table 
[2] The correct frequency is marked with a dashed vertical line. The highest frequency (/ ~ 0.136J -1 ) corresponds to a period just 
twice that of the binary - near the instability limit. The periodograms are normalized to the their own RMS (as described in Paper 
I) - thus, a value of 5 is a 5<x detection of a certain period. 



of constant stars converges. Note that in each iteration the effects 
must be always computed from the same raw data and not from 
the SYSREM-corrected data. This procedure also improves the 
signal-to-noise ratio of the SYSREM effects. For the TrES Lyrl 
data we used the final set of constant stars to calculate the three 
SYSREM effects, with ultimate precision similar to that of the 
TrES survey. 

2. Detrending variable stars: The constant stars allowed us 
to determine the systematic effects well - and from now on the 
effects are assumed known and will not be changed. To correct 
variable stars for these systematic effects we will have to apply 
the same correction determined above to the residuals around the 
smoothed LC of the target variable star by iteratively smooth- 
ing the target LC, detrending the residuals, adding back the cor- 
rected residuals to the previous smooth and then re-smoothing 
the LC, till convergence While this procedure is sensitive to the 
exact smoothing technique, when the smoothing parameters fit 
well for a given LC the results are quite significant for transit 
detection (see Fig. [7]), and allow us to reduce the scatter of the 
LC around the model to the same level as the scatter of LCs of 
constant stars of similar brightness. 

3. EB solution: Light curves for eclipsing binaries in the 
Lyrl field were analyzed by the NDE model. Errant observations 
were picked out by eye and eliminated before analysis. Accurate 
dates on minima were determined from the data by using the 
method of Kwee & van Woerden (1956). The resulting periods 



and zero epochs are listed in Table [3] along with the fitted orbital 
parameters. In most cases, the fitted light curves were relatively 
insensitive to the assumed values of the limb-darkening coeffi- 
cients (x p and x s ) and to the mass ratio (q), which was guessed 
based on preliminary values of the light ratio. 5 of the 17 light 
curves showed significant eccentricity. 

4. Generation of final EB LC: Almost identically to the 
"Detrending variable stars" step - the raw LC is detrended us- 
ing the solved EB model. Now the model is final so multiple 
iterations are not needed. 

Fig. [8] depicts further examples of EBs from the TrES Lyrl 
field, showing that the generation of multiple near-systematics- 
free EB light curves from wide-field transit surveys using the 
procedure above is indeed possible. The low level of systematics 
will help CB-BLS to be less likely to find spurious signals. All in 
all, about two dozen detached or semi-detached EBs were bright 
enough (similarly-bright constant stars had ultimate precision of 
~ 2%) to be processed as above. Of the two dozen EBs iden- 
tified 17 were fitted, and for 8 of them planet detection seemed 
possible (the other EBs exhibit either low inclination angles, evi- 
dence of spots or period change), and CB-BLS was applied - but 
we did not detect any significant signal. The acB = 12 detection 
limits for each of the objects are given in Fig. [9] 

We comment that since there were about 7000 objects in the 
relevant magnitude range one may expect about 70 EBs to be 
present in the data, and thus a yield of two dozen EBs may seem 
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Fig. 5. The Test3 LC residuals folded against scaled-distance for the best-fit model (see text for details). Readers who are familiar 
with phase-folded LCs can view of this figure as yet another fold of the LC about the center of the transit, so that both ingress and 
egress fall on scaled-distance ~ 1 , and all in-transit points are expected below scaled-distance of 1 (to the left of the vertical dashed 
line). To aid visibility the data is binned to 0.1 units of this scaled-distance. Despite the fact that Test3 was the shallowest signal in 
this blind tests series CB-BLS clearly separates in-transit from out-of-transit points. 



Table 3. Results of fitting light curves of bright EBs from the TrES Lyrl field. Boldfaced star names (NStED ID) indicate that 
CB-BLS was subsequently applied. The other columns are: the period P and reference epoch E Q (in HJD-2453000), the surface 
brightness ration J s , the primary, secondary radii and their ratio r p , r s and k, the orbital eccentricity, inclination and argument of 
periastron passage i, e and oj, the light fraction of third light, primary and secondary components L3, L p and L s , the linear limb- 
darkening coefficients x p = x s , the mass ratio q = the RMS of the data around the model cr, and the number of data points used 
in the fit N. 
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too low. However, a similar number of EBs were actually iden- 
tified but prioritized very low and not processed further (such 
as EBs with long periods, grazing EBs, EBs with low SNR, 
EBs with additional signals overlaid, non-eclipsing contact EBs, 
etc.). Considering that further EB's were missed due to the win- 
dow function - the total fraction of EBs seems roughly compati- 
ble with CoRoT's ~ 1% above. 



4.2. Other targets 

While wide-field transit surveys such as C0R0T will probably be 
the main source of EBs for future CB-BLS analyses (as demon- 
strated above), there are currently a number of open literature pa- 
pers on EBs that have a long enough and accurate enough LCs 
to possibly allow the detection of transiting planets. These pa- 
pers focus on the EBs and so have the added benefit of having 



the EBs (usually) fully analyzed - including spectroscopic de- 
termination of the mass ratio q (which allows us to remove one 
free parameter, and shorten the CB-BLS run time). Note that 
for some of these objects the new multi-band CB-BLS feature 
is required. Also note that while the blind test were conducted 
on a ~ 50 d long LC, these objects were observed over much 
longer periods - typically several years. This necessitates a cor- 
respondingly higher frequency resolution, and consequently - 
significantly longer CB-BLS run time. 

Light curves of WW Cam (Lacy et al 2002) and V1061 Cyg 
(Torres et al. 2006) were taken from the literature and fitted by 
the NDE model (Etzel 1981, Popper & Etzel 1981) as described 
above. LC parameters for those stars are listed in the referenced 
articles. Light curves of V432 Aur [Siviero et al. 2004] and BP 
Vul [Lacy et al. 2003] were also taken from the literature and fit- 
ted with JKTEBOP - giving LC parameters similar to the ones 
listed in the referenced articles. In FigfTOlwe plot the above LCs 
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Fig. 6. The Alarm statistic for 1000 objects (all objects in some 
magnitude range) of the TrES Lyrl dataset. The insert shows 
the lower part of the Alarm distribution before the first iteration 
for variable-stars filtering (see text). After we used only stars 
with Alarm< 35 to calculate one SYSREM effect - the Alarm 
distribution of the same objects looked as shown on the main 
panel. 

TrES Lyr1 02432 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



Fig. 7. Two light curves of star 02432 from the TrES Lyrl field. 
Bottom, Red LC: the light curve as posted on the public NStED 
database: it was detrended with no special care for intrinsic vari- 
ability. Top, Blue LC: the same raw data processed as described 
in the text - with cleaning procedures applied only to the resid- 
uals. Note that transit detection in LCs such as the red LC is 
nearly impossible. 



together with their best-fitting model and model residuals. We 
applied CB-BLS to these systems too but did not detect any sig- 
nificant signal. 



5. Conclusions 

We presented a number of improvements to the CB-BLS algo- 
rithm that allow CB-BLS to improve it's sensitivity to shallow 
transits, correctly model distorted close binary stars, use multi- 
band photometric data, and more. We also presented an analyt- 
ical approach for the determination of detection limits for CB- 
BLS. The resultant implementation of CB-BLS was blind-tested 
using five test light curves derived from CoRoT data — all five 
were correctly and strongly identified, and the theoretical pre- 



Fig.8. More examples of the EBs from the TrES Lyrl field. 
These EB LCs were generated using the process discussed in 
the text for the proper detrending of variable stars. This, together 
with good modeling, resulted in extremely small systematics in 
the residuals across the > 11, 000 data points. The object ID at 
the top of each panel is the NStED designation. 
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Fig. 9. The acB = 12 detection limits for each of the TrES EBs 
marked in Table [3] — i.e., the smallest depth of primary transits 
caused by a hypothetical CB planet in the systems that is ex- 
pected to be reliably detected w ith C B-BLS. The detection limit 
was calculated as explained in §2.2| The planetary periods used 
are from twice the host EB period to half the survey time span. 



diction for the detection limit was in line with the actual results. 
We also presented an example of a procedure for the correct de- 
trending of variable stars from wide-field surveys that results in 
far less systematic input light curves for CB-BLS, and a subse- 
quent first application of CB-BLS to real LCs from both targeted 
observations and a wide-field survey. 

Looking for transiting CB planets requires good understand- 
ing of several issues: how to generate good photometry of 
intrinsically-constant stars (now well understood), how to gen- 
erate good photometry of intrinsically- variable stars (for e.g., 
as explained above: iterative detrending of the residuals), high- 
quality EB modeling, and finally - the actual search for transiting 
CB planets (e.g., the CB-BLS algorithm). While the above gen- 
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Fig. 10. Some of the open literature EBs to which the CB-BLS 
analysis was applied. Each EB light curve is phased to it's or- 
bital period with the best-fit model over plotted and the model 
residuals drawn below. 



eral description will probably remain correct, we are still in the 
process of learning the details of the process - as exemplified by 
the still-evolving CB-BLS. Still, we feel that the experience of 
the single-stars transit surveys allows for an even steeper learn- 
ing curve (and we remind the reader that single- stars transit sur- 
veys had their own learning curve: more than 80% of the known 
transiting planets today were not known only three years ago). 

Looking for CB planets around EBs is advantageous since 
EBs are already well-aligned to our line of sight. Since theo- 
retical models predict that CB planets will be, at least initially, 
aligned with their host EB - the chances that CB planets will 
transit in front of the EB components are relatively high. CB 
planets in general and transiting CB planets in particular are ex- 
pected to be particularly beneficial for the general study of ex- 
trasolar planets: 

- Binaries are a very significant fraction of the total stellar pop- 
ulation. Studying planet formation and evolution without in- 
cluding the formation and evolution of planets in binary star 
systems is surely incomplete. 

- Frequency of CB planets depends on assumed planetary for- 
mation mechanism, and may allow to distinguish between 
competing theories. The detection limits formalism intro- 
duced here will allow us to give meaningful constrains on 
the frequency of CB planets by correctly accounting for null 
results. 

- After the initial surprises of close-in giant planets and plan- 
ets in eccentric orbits, detecting CB planets will again sub- 
stantially diversify the environments that produce and sus- 
tain planets. 

- Understanding very high precision LCs of transiting plan- 
ets around single stars is limited by the uncertainty of stel- 
lar parameters [e.g., Johnson et al. 2008, Gillon et al. 2008, 
and discussion therein]. EBs allow to us significantly in- 
crease stellar parameters accuracy, so resultant fits to plan- 
etary structure models will have smaller error bars. 

- Follow up will be interesting for several reasons: I) close-in 
CB planets can show orbital evolution on short time scales 
(even 100s of days), II) CB planet systems produce four dis- 
tinct Rossiter-McLaughlin [e.g., Gaudi & Winn 2007] ef- 
fects: 1-2, 2-1, p-1 and p-2 (where 1,2 and p designate the 
two binary components and the CB planet respectively, and, 



for e.g., 2-1 stand for the spectroscopic eclipse of component 
1 by component 2), III) Higher chances of finding resonant 
systems and/or chaotic systems. 

Specifically, using CB-BLS to look for transiting CB planets 
is especially useful since: 

- CB-BLS allows to find shallow transiting CB planets in the 
residuals of EBs. 

- CB-BLS can naturally include Roche geometry - and there 
are many EBs that require that. 

- CB-BLS can simultaneously fit multi-band data. 

- CB-BLS allows to harness existing datasets to the detection 
of transiting CB planets with no further inputs. Medium (or 
maybe even low) resolution spectroscopic input can allow to 
eliminate the mass ratio as free parameters and to reduce the 
required CPU time, but this is not mandatory. 

- CB-BLS source code and documentation are planned to be 
made freely available. The source code is almost entirely in 
MATLAB and so easily portable. 
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